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Abstract 

We study the behavior of two-dimensional Dirac fermions in the presence of a static 
long-range-correlated random vector potential. By applying an exact path integral 
representation for the propagator of a spinor particle we obtain asymptotics of the 
gauge invariant spectral function and the correlation function of the local density 
of states, both in the ballistic regime of sufficiently high energies. We also discuss 
localization properties of the random Dirac wave functions in the complementary 
zero energy limit and the putative localization scenario. 
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1 Introduction 



In recent years, the standard theory of electron localization in two dimensions 
has been extended to the situations where, instead of potential disorder, the 
fermions are subject to a random vector potential (often referred to as random 
magnetic field or RMF) whose Gaussian spatial fluctuations are described by 
a pairwise correlator 

(a i (q)a j (q')} =w ij (q)5(q + q), (1) 
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where w^q) = ((% - 

This problem was first encountered and then extensively studied in the con- 
text of the compressible Quantum Hall Effect at even denominator filling frac- 
tions. In order to describe the ostensibly Fermi liquid-like properties of these 
strongly correlated electronic states, a new kind of quasiparticles ( "composite 
fermions") which would effectively "see" an ordinary potential disorder as a 
long-range-correlated RMF described by Eq.(l) with 77 = 2 was introduced in 
[1,2]. 

Shortly thereafter, a (pseudo)relativistic version of the short-range-correlated, 
77 = 0, RMF problem had been proposed as the means of describing the Integer 
Quantum Hall plateau transitions [3] where it was later identified with the 
continuum limit of some network models. The intrinsic scale invariance of the 
77 = problem prompted the use of the powerful machinery of two-dimensional 
conformal field theory in most of the subsequent work on this topic. 

Recently, a novel (this time, long-range correlated, 77 = 2) variant of the rela- 
tivists RMF problem has emerged in the theory of localization of the Dirac- 
like nodal quasiparticles in the mixed state of planar <i-wave superconductors. 

In the presence of vortices of the <i-wave superconducting order parameter, not 
only do the nodal quasiparticles' energies get semiclassically Doppler shifted 
due to the circulating supercurrent [4], but also their wave functions acquire 
intrinsically quantum mechanical Bohm-Aharonov (BA) phases. 

In order to account for both effects on equal footing, a singular gauge trans- 
formation similar to that of Ref.[2] was implemented in Ref.[5]. In the case 
of a disordered vortex array, which corresponds to the experimentally well- 
documented vortex line liquid phase where the vortices are believed to be 
randomly pinned by columnar or other strong defects, the transformation of 
Ref. [5] converts the nodal quasiparticles into auxiliary neutral Dirac fermions. 
The latter, in addition to a random scalar potential accounting for the Doppler 
shift, are subject to an effective RMF with zero mean described by Eq.(l) 
where 77 = 2 and a is proportional to the areal density of vortices [6,7,8,9]. 

Moreover, a formally similar RMF problem has emerged in the theory of the 
pseudogap phase of the cuprates described as a plasma of thermally excited 
vortex-antivortex pairs, the parameter a being proportional to temperature 
[10,11]. Furthermore, one encounters yet another formal analog of this problem 
when analyzing the effect of topological structural defects (dislocations) on the 
Dirac-like electronic excitations with the momenta near conical ( Appoints in 
the hexagonal Brillouin zone of a graphite layer [12]. 

However, despite its occurrence in all of the above physical situations, the 
77 = 2 RMF problem for massless Dirac fermions has, thus far, remained 
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unsolved. In fact, this problem does not appear to be readily amenable to any 
of the methods applied to its previously studied rj = 2 non-relativistic [1,2] 
and rj = Dirac [3] counterparts. 



In the present paper, we fill in this gap by developing a novel non-perturbative 
path integral approach to the single-fermion action 



S = J dr J dtip(t,r) i%d t — 7*(V — a(r)) ip(t,r 



(2) 



where the correlation properties of the static random vector potential a(f) 
are governed by Eq.(l), and the 7^-matrices obey the algebra {7^,7^} = IS^ 
(hereafter, we put the speed of the Dirac fermions equal to unity). 



2 Path integral representation for gauge-invariant propagator 



In order to gain a preliminary insight into the problem we first attempt to 
apply the customary self-consistent Born approximation (SCBA) to the con- 
ventional (retarded) fermion propagator G R (t, r) = Q(t){ip(t, r)ip{0, 0)). A 
closed SCBA equation can be obtained for the fermion self-energy which is 
defined via the Fourier transform G R (e,p) = [e% — pj + S R (e,p)] _1 

Z R (e,p) = J T^mMliG^P + q)^- (3) 



Although one can readily solve Eq.(3) for any 77 < 2, including the previously 
studied case of rj — [3], at rj — 2 the analysis is impeded by the fact that this 
equation now exhibits an infrared divergence with the size L of the system 
(E oc 7oVln L), thus rendering the SCBA inapplicable. 

This intrinsic divergence remains for all e and p, even if one proceeds beyond 
the SCBA by replacing (3) with the equation E = nT*G, where n stands for 
the density of Poissonian-distributed localized RMF sources (hereafter, moti- 
vated by the discussion of the mixed state of d-w&ve superconductors, we refer 
to them as "vortices"). Although the exact T-matrix given by the expression 
T = w * (I — G* wy 1 does account for all the events of multiple scattering 
by the same vortex, the above infrared divergence can not be made disappear, 
for it is directly related to the non-gauge invariant nature of G R (t,r). 

This problem does not arise, however, when calculating such observables such 
as the gauge-invariant Green function of the nodal quasiparticles introduced 
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in in the context of the d-wave superconductors [13,14,?] 

Gt(t, r) = B(t)m, r) exp ^-i J 3(r*)dr*j ?(0, 0)), (4) 

which is manifestly invariant under arbitrary time-independent gauge trans- 
formations for any choice of the contour T. 

In Refs. [13,14,15], where the amplitude (4) was discussed in the case of a 
dynamic vector potential a(t, f) (the latter being a part of the fictitious time- 
dependent gauge field representing either spin [13] or pairing [14,15] fluctu- 
ations), the contour T was chosen as a straight line between the end points 
(0,0) and (t,r). 

It must be noted, however, that in the Lorentz-invariant situation considered 
in Refs. [13,14,15] the amplitude (4) with the straight-line contour T appears to 
exhibit an unphysical power-law behavior characterized by a negative (instead 
of a positive, as one would have expected on the general physical grounds and 
as Refs. [13,14,15] erroneously claimed to be the case) anomalous exponent 
[20]. 

Although, thus far, no physically sound alternative to the heuristic form of 
Eq.(4) conjectured in Refs. [13,14,15] has been found, it turns out that in the 
non-Lorentz-invariant case of a static vector potential a(r) the situation is 
rather different. 

Namely, in the case of interest, the most important features, such as the mean 
free path (see Eqs. (11,12) below) which controls the exponential (rather than 
a power-law, as in the case of a dynamic RMF [13,14,15,20]) decay of Eq.(4) 
or the structure of the "near-shell" singularity of the fermion spectral function 
(see Eqs. (13,14) below), remain robust against deformations of the contour T 
(apparently, in the static case the only potentially important is the spatial 
projection of T, whereas its pitch in the temporal direction is not). 

Moreover, while being representative of the properties of a generic gauge- 
invariant one-fermion amplitude, the above choice of the contour T facilitates 
the calculation of Eq.(4) along the lines of the earlier analyses of the non- 
relativistic version of the rj = 2 RMF problem [16,17,18,19]. 

Of course, the spinor nature of the amplitude (4) makes it impossible to apply 
the results of Refs. [16,17,18,19] directly. Therefore, we first utilize the (rela- 
tively unknown) path integral representation devised in [21] and, for a fixed 
RMF configuration {a(r)}, cast the time Fourier transform of Eq.(4) in the 
form of a path integral over the spatial coordinate r (r) and an auxiliary vector- 
like variable k{r) parametrized by the proper time r (not to be confused with 
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the real time t\) 



oo r(r)=r 

dr 

r(0)=0 



G;';, v (f. r\a) > / dr / DrD/,- exp | /S (Y) + i / a(r)^dr'-i / ^(f')df' ] (5) 



where the first term in the exponent represents the (matrix- valued) free fermion 
action 



So(r) 



dr' 



(6) 



while the other two correspond to the RMF coupling term in Eq.(2) and the 
line integral in the exponential factor inserted into Eq.(4), respectively. In 
Eq.(5), the usual r-ordering of the 7-matrices must be performed according 
to the order of their appearance in the series expansion of the exponent. In 
contrast to various approximate (e.g., Bloch-Nordsieck) representations, the 
integration over k{r) allows one to exactly account for the spinor structure of 
the fermion propagator. 

Averaging (5) over the Gaussian RMF introduces the exponential attenuation 
factor 



W[f(r)] = exp 



AT) 



whose argument, in the case of 77 = 2, turns out to be proportional to the 
so-called Amperian area of a closed contour which is composed of a given 
trajectory r(r) and the straight-line segment traversed backwards (from r to 
0) [16,17]. 



3 Gauge-invariant propagator in the ballistic regime 



While being strictly positive, the Amperian area coincides with the absolute 
value of the regular (algebraic) one only for non-self-intersecting closed con- 
tours, which, strictly speaking, precludes Eq.(7) from being interpreted as 
an additional term in the effective Dirac action. Nevertheless, this becomes 
possible for the trajectories whose projections onto the straight path Y are 
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single-valued, whence (7) reduces to 

W[r(r)}^exp(-^j \x ± \dxA , (8) 



where the parallel and transverse components of the vector r(r) = (x\\,x±) 
are defined with respect to the direction of the fermion's propagation. 

The relativistic kinematics of the Dirac fermions makes the present situa- 
tion markedly different from the previously studied non-relativistic version of 
the f] = 2 RMF problem [16,17,18,19]. As opposed to the interpretation of 
the non-relativistic analog of SS = i\nW[f(r)] as an effective linear confin- 
ing potential [16,17,18,19], Eq.(8) can be best thought of as an (imaginary) 
position-dependent mass whose linear increase with restrains the Dirac 
fermion's motion in the direction perpendicular to its "classical trajectory" 
r (r) = r(r/r), the latter representing the saddle point of the path integral 
(5). 

Comparing the saddle-point action S ~ er to the characteristic RMF correc- 
tion 5S ~ (ae/pY^r, we find that the condition So ^> 5S, under which the 
path integral (5) is dominated by the trajectories close to ro(r), is satisfied in 
the "ballistic" regime of high energies and large momenta (ep ^> a) . 

In this regime, one can integrate over the momentum p(r), thereby effectively 
squaring the Dirac operator (for details, see [20,21]), as well as the longitudinal 
spatial x\\(t) component of the coordinate vector r(r), thus reducing (5) to 
the one-dimensional path integral over x±(t). 

As in the non-relativistic rj = 2 RMF problem [16,17,18,19], the computation 
of the remaining path integral amounts to finding the resolvent of the second 
order ordinary differential equation 



{pl - e 2 ) - d]_ + (ae/pA x 2 ± + ia sign x± 



g(e,p\\\x±,x' ± ) = 5(x± - x' ± 



whose exact form can be readily obtained in terms of the Wronskian of a pair 
of linearly independent solutions of the corresponding homogeneous equation 



#(e,P|||0,0) 



E 



r(z A + 1/ 2) 
r(*A) 



-i -i 



(10) 



where z\ = ip\\(e — \p\\)(e + \p\\ + ia/p\\)/Aea and T(z) is the standard T- 
function. 
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The asymptotical large-distance behavior of the gauge-invariant propagator 
is controlled by the closest to the real pn-axis zero of the denominator in 
Eq.(lO). Applying to g(e,p\\\0,0) the one-dimensional Fourier transformation 
with respect to p\\ we arrive at the expression 

G£ v (e,f) ~ ( e7o - —7] -^==sin(|e|r + 0i), (11) 



where 0i ~ 1, which decays exponentially at r — > 00 with the characteristic 
length scale 

/ w ^1 (12) 
71a 



The latter can be naturally identified with the RMF mean free path which, 
according to Eq.(12), decreases as the fermion energy gets lower, contrary 
to the situation for rj = [3] (but qualitatively similar to that found in the 
non-relativistic variant of the rj = 2 problem where / ~ lel^a-Va [16,17,19]). 

The "near-maximum" (|e 2 — p 2 \ < a) behavior of the electron spectral function 
is determined by the two-dimensional Fourier transform of ImG? v (e, r) which 
can be obtained by convoluting Eq.(10) with the kernel (p 2 — g 2 )~ 3//2 , thus 
resulting in the expression 

lmGt(e,p) « V^f( e 2_I?2/J/n2i3/4 » ( 13 ) 



which explicitly demonstrates that the bare pole gets replaced by branch cuts 
of the function z^/ 2 (hence the emergence of the exponent 3/4 in (13)). 

At larger deviations from the maximum (a <^ |e 2 — p 2 \ < e 2 ) the spectral func- 
tion is determined by the large- z\ asymptotics of the T-functions in Eq.(10). 
In this limit, the spectral function resembles the Lorentzian 

ImG inv (e,p) oc a— 14 

[e A — p z Y + a z 



although, under a closer inspection, the shape of ImG^ v (e, p) turns out to be 
asymmetrical as a result of the Lorentz invariance's being broken by the static 
disorder. Moreover, the product of its maximum height and width reveals a 
marked departure from the Lorentzian value (~ 7r/4 instead of 1). 
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4 Averaged products of Green functions 



Our eikonal-type approach also enables one to study other gauge-invariant 
physical observables, such as a correlation function of the local density of 
states (DOS) which is related to the average (G R G R ). Owing to its gauge 
invariant nature, the latter is free of infrared divergencies and can be cast in 
the form of a two-particle path integral 



oo oo r r 

(G R (e,r)G R (e, -f)> = J dn J dr 2 J Dr&h J Dr 2 Dk 2 



e ,S (n) e ,So(r 2 ) -Q Wifiin)-^)), (15) 
i,J=l,2 

where the product of the lU-factors yields the exponent of the Amperian area 
of the closed contour formed by a pair of trajectories fi(r) and r 2 (r). 

We mention, in passing, that the average (15) equals that composed of a pair 
of gauge-invariant amplitudes ({G R G R ) = (G R V G R V )), since, regardless of the 
choice of T in Eq.(4), the phase factors from the two functions G^ v exactly 
cancel against each other. 

In the ballistic regime, the integral (15) receives its main contribution from the 
pairs of trajectories with single- valued projections to the direction of the vector 
f. Upon separating the coordinate variables onto the "center of mass" and the 
relative motion (r* = r\ ±r* 2 , = h\ ±/c 2 ) parametrized by = T\ ±r 2 and 
integrating over all the variables but xj(r + ), one finds that the average (15) 
can be again related to the solution of Eq.(9), albeit with some extra powers 
of two stemming from the Jacobian of the above transformation. As a result, 
the variance of the local DOS (5v(e, r) = z/(e, f) — {v{e))) decays with distance 
as 

(5v(e, r)6u(e, 0)) oc 1 1 sin(2|e|r + 2 ) , (16) 

y/rl 



where 2 ~ 1. 

Notably, the average of the product of two Green functions (16) does not 
amount to the product of the two averages (each of which is given by Eq.(ll)), 
thereby indicating the presence of non-trivial vertex corrections. 

Nevertheless, the exponential decay of Eq.(16) is controled by a length scale 
equal to the half of the mean free path (12), which corroborates our conclusion 
that (alongside a number of other robust features such as the "near-shell" 
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asymptotic behavior (13)) Eq.(12) is, in fact, independent of the choice of the 
contour T in Eq.(4). 

In a similar manner, one finds the average {G A G R ) = (G^ V G^ V ) oc (e 2 /7r) 1 / 2 e~ 2r / / 
which can be used for computing such transport coefficients as optical con- 
ductivity a(u,T). However, such a result would be limited to the regime of 
high frequencies or temperatures max(w,T) ^> a 1 / 2 where other mechanisms 
of scattering, which may exist in the known realizations of the RMF problem 
(see Refs. [6,7,8,9,10,11,12]), might become important. 



5 Localization of zero-energy modes 



At low energies and small momenta ep < a our approach ceases to be valid, 
and, for instance, the low-energy behavior of the DOS can not be readily 
inferred from the above discussion. In order to shed some light on the effect of 
RMF scattering on the low-energy part of the Dirac spectrum we look into the 
properties of the zero-energy states. Quite remarkably, the latter exist for an 
arbitrary configuration of the non-uniform magnetic field, thanks to the exact 
cancellation between the orbital and the Zeeman terms in the total energy of 
a spinor particle with the gyromagnetic ratio equal two. 

Parametrizing an arbitrary (up to a gauge transformation) RMF configura- 
tion in terms of a scalar function Oj(r) = €ijdj&(r) one can write a pair of 
independent zero-energy wave functions (there are no other states with e = 
as long as the total RMF flux vanishes) as follows 

/ p*(«0 \ 

*±(rO«(l±7o)( e _» (f) J • (17) 



In a finite-size system, the degree of the wave functions' localization (or a lack 
thereof) can be inferred from the set of inverse participation ratios 

Jl*(r)l 2 Mf 

n ^L 2 (/|*(r*)| 2 dT s ) n/ ' 1 ' 

Performing the Gaussian average over the disorder field $(r) with the weight 
P[$(f)] oc exp(- / df(V 2 $) 2 /2a:) we obtain 

Pn OC — , (19) 

which is suggestive of strong localization of the zero-energy states, the "local- 
ization length" being of order £ ~ a -1 / 2 . This conclusion is in stark contrast 
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with the i] = case [3] where the prelocalized zero-energy wave functions ex- 
hibit a multifractal spectrum of anomalous dimensions P n oc L~ an+bn . Con- 
ceivably, as the energy decreases past a 1 / 2 , the elastic mean free path (given 
by Eq.(12) at high energies) saturates at I ~ £, followed by the onset of local- 
ization at still lower energies. 

In the f] = case, it was recently argued [22,23] that the corresponding lo- 
calization scenario belongs to the so-called "C" universality class which en- 
compasses the ci-wave superconductors with a strong inter-node scattering 
in the absence of time reversal symmetry [24,25]. It was also conjectured in 
Refs. [22,23] that the low-energy properties of the system are described by the 
non-linear a-model (NLcrM) on the coset space OSp(2|2)/GL(l|l). 

In order to extend this conjecture to the Dirac rj = 2 RMF problem one would 
have to derive an effective NLcrM of the appropriate symmetry by carefully 
separating between the massless and massive modes of the (this time, non- 
local) Hubbard-Stratonovich disorder fields, by analogy with the procedure 
carried out in the nonrelativistic 77 = 2 RMF problem [26]. 

There is a good reason to believe that the resulting localization scenario falls 
into the already existing classification of the ci-wave universality classes char- 
tered in Refs. [24,25]. However, given the notorious sensitivity of the ci-wave 
localization patterns not only to the type (potential vs magnetic vs extended 
defects, etc.) but even to the strength (Gaussian vs unitary) of disorder, such 
important details as the asymptotical form of the low-energy DOS ((^(e)) oc e 2 
for a generic class "C" system [27]) or a possible crossover from "C" to an- 
other, e.g., the so-called "A", universality class at intermediate energies due 
to the predominantly small-angle nature of the BA scattering, remain to be 
worked out. 



6 Summary 

To conclude, in the present paper we investigated the properties of the two- 
dimensional Dirac fermions subject to a long-range-correlated random vector 
potential. In the ballistic regime of large quasiparticle energies, we obtained the 
asymptotics of a representative gauge-invariant fermion Green function and 
the DOS correlation function. In the complementary low-energy regime, we 
found a signature of strong localization of the zero-energy states by computing 
their inverse participation ratios. 

Thus far, it has been rather difficult to find a conclusive experimental evidence 
of the quasiparticle localization in the high-T c cuprates. However, theoretical 
predictions based on the results of this paper which pertain to the more exper- 
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imentally accessible ballistic regime may be possible to test with a number of 
standard probes (thermal transport, ARPES, tunneling, specific heat, NMR) 
performed in the vortex line liquid phase of the superconducting cuprates [28]. 

To this end, it is worth mentioning that the results of a recent numerical 
analyis of the low-field thermal quasiparticle conductivity in the vortex line 
liquid state of the ci-wave superconductors [29] appear to be in good agreement 
with the analytical result based on the analog of Eq.(12) [28]. 

Moreover, the results of this paper can be used to study the effects of both 
thermal [10,11] and quantum [13,14,15] phase fluctuations on the spectrum of 
the nodal quasiparticles in the pseudogap phase of the cuprates as well as that 
of randomly distributed dislocations on the Dirac-like electronic excitations in 
graphite [12]. 

In particular, when applied to the theories of the pseudogap phase [13,14,15] 
where the nodal fermions are scattered by the quasi-static (thermal) gauge 
fluctuations which represent either the effect of a disordered flux-phase [13] or 
that of thermally excited vortices [14,15], the results of this paper imply that 
the gauge-invariant propagator should exhibit an exponential, rather than a 
power-law, decay, contrary to the conclusions drawn in Refs. [13,14,15]. 
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